Numerical Evolution in time of curvature perturbations in Kerr black holes 
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In this paper I will review the basic features of the theory of curvature perturbations in Kerr 
spacetime, which is customarily written in terms of gauge invariant components of the Weyl tensor 
which satisfy a perturbation equation known as the Teukolsky equation. I will describe how to evolve 
generic perturbations about the Kerr metric and the separable form of the wave solutions that one 
obtains, and the relation of the Teukolsky function to the energy of gravitational waves emitted by 
the black hole. A discussion of a numerical scheme to evolve perturbations as a function of time 
and some preliminary results of our research project implementing it for matter sources falling into 
the black hole is included. 



I. INTRODUCTION 



The advent of the possibility of detecting the existence of gravitational waves with large interferometers (like 
the LIGO or VIRGO projects now in their final construction stages) is a very exciting prospect for relativists and 
astrophysicists. The weakness of the predicted signals requires the search for very strong sources. The best candidates 
to date for detectable sources are black hole-black hole collisions. There is a very intense effort to numerically solve 
Einstein's equations in full non-linear form to provide accurate templates to aid in the detection of the waveform 
from such a collision. Due to the complexity of this task, this is a very ambitious problem still untractable even by 
' state-of-the-art supercomputers. Approximations like treating the last stages of the collision as a perturbation from a 
single stationary black hole spacetime are very helpful both for aiding in the construction of a code for a full non-linear 
■ evolution and to shed light on the physics of the collision itself. 

In black hole perturbation theory we want to linearize the Einstein equations for gravitational (and other matter 
fields) about one of the known stationary solutions. This is a simplification over solving the full non-linear set of 
', equations, although the resulting equations can be quite messy in algebraic terms. But in this way we can evolve small 
departures from the isolated black hole spacetime due to the presence of outside sources like particles, gravitational 
qh| waves or even collisions with other black holes (in the "close limit" where the two objects can be regarded as one 
^ 1 distorted black hole). 

A straightforward way to do this is to consider metric perturbations. This is analogous to what is done in the 
linearized treatment of gravity when the metric is written as g a b = r) a b + h a b, where h a b is small compared to the 
. *-h \ Minkowski metric. In the late 50's, Regge and Wheeler [Q were the first to do this by linearizing the vacuum Einstein 
equations about the Schwarzschild metric and using the symmetries inherent in the background metric to separate 
the angular and radial parts of the equation. They expand the solution in tensor spherical harmonics, leaving a rather 
simple Schrodinger-like equation to solve for the radial part of the odd-parity set of perturbations. 

The Regge- Wheeler approach, even when it has been quite useful for the Schwarzschild case, involves quite a 
lot of complex algebra, particularly when treating even-parity perturbations. It is unfeasible to extend it to the 
more complex Kerr geometry. In addition, it suffers some of the same drawbacks as the linearized theory, since the 
perturbations are coordinate gauge dependent. If one does an infinitesimal coordinate translation x a — > x a + e£ a , the 
metric function h a b transforms as 

hab -> Kb - 2V (b £ a) (1) 

The use of gauge dependent perturbations is warranted if one asks the right questions about global properties of the 
gravitational radiation and avoids trying to say anything about properties of the gravitational field at local events in 
the spacetime. Still one would be able to avoid all ambiguities if only gauge independent perturbations are considered. 



A. Special algebraic character of the Kerr metric 

The Kerr metric [Q] is the axisymmetric solution to Einstein's field equations that describes the spacetime outside a 
stationary, rotating black hole. It is a type II-II in the Petrov classification scheme, which means that out of the four 
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principal null directions that all stationary spacetimes can have it has two distinct principal null directions, (since the 
other two coincide with these). 

These directions are null vectors that satisfy the following condition : 

k b k c C abc[d k e] = (2) 

where the C abcd is the Weyl tensor. 

This fact will be very convenient when one treats curvature perturbations using the Newman-Penrose formalism 
||, since if one selects these principal null directions as the basis for the NP tetrad, one gets enough relations between 
the various spin coefficients to make the resulting system of equations easily solvable. 



B. The Newman- Penrose formalism 



The Newman-Penrose formalism was developed to introduce spinor calculus into general relativity. It is a special 
instance of tetrad calculus Q] . Let us present a brief summary of the basic ideas behind the use of the formalism 
which will be used in deriving the Teukolsky equation. 

One starts by introducing a complex null tetrad { I, n, m, m* } at each point in spacetime which consists of two 
real null vectors I, n and one complex spacelike vector m. These should satisfy the orthonormality relations 

I ■ n — m ■ m* = 1 (3) 

with all other products being zero. Then g ab = — 2[/( a nt) — m( m£\]. Now one defines four directional derivative 
operators along the tetrad directions 

D = l a V a ,A = n a V a , 

6 = m a W a ,5*=m* a W a . (4) 

The basic quantities of the formalism are the spin coefficients, of which there are twelve complex ones : 

a = | (n a m* b X7 b l a - m* a m* b \7 b m a ), 
p= \ (n a m b W b l a - m* a m b W b m a ), 
7=| {n a n b V b l a - m* a n b W b m a ), 
e=\ (n a l b V b l a -m a l b V b m a ), 
A = —m* a m* b 'V b n a , /j, — —m* a m b 'V b n a , 
v = -m* a n b V b n a , n = -m* a l b 'V b n a , 
K = m a l b V b l a , p = m a m* b V b l a , 

a = m a m b V b l a , t = m a n b V b l a (5) 

The whole set of field equations in the formalism come by writing the Ricci and Bianchi identities using these 
coefficients, and they take the place of the Einstein equations. 

All ten independent components of the Weyl tensor can be written as five complex scalars: 

<Ao = -C abcd l a m b l c m d 
<Ai - -C abcd l a n b l c m d 
%l) 2 = -\ C abcd {l a m b l c m d + l a n b m c m* d ) 
i>3 = -C abcd l a n b m* c n d 
i>A = -C abcd n a m* b n c m* d (6) 

To do perturbation calculations one specifies the perturbed geometry by introducing slight changes in the tetrad 
like I — l A + l B ,n — n A + n B , etc. Here the A terms are the unperturbed values and the B ones the small 
perturbation. Then, all the Newman-Penrose spin coefficients and other quantities can be also written in a similar 
fashion : ip4 — i\) A + ip^ , etc. The perturbation equations come from the Newman-Penrose set by keeping B terms 
only up to first order. 
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II. THE TEUKOLSKY EQUATION 

In 1973, Teukolsky || used the Newman-Penrose formalism to the special case of the background geometry of a 
II-II type, (the Kerr or Schwarzschild black holes are of this type) . In this way he was able to deduce the linearized 
equations for full dynamical perturbations of the hole that could handle changes in its mass and angular momentum, 
interaction with accreting test matter or distant massive objects, etc. 

This approach has many important advantages. First, it turns out rather surprisingly that the equations are 
separable, so that by Fourier transforming and expressing the solution as a series expansion of spheroidal harmonics 
one ends up with having to solve just an ordinary differential equation for the radial part just like in the Regge- Wheeler 
case (in fact, the solutions are related to each other by a transformation operator, as we will discuss later). Second, 
since for gravitational perturbations the dependent variable will be constructed out of the Weyl tetrad components 
■00 an d ipA j this will describe gauge independent perturbations, because these are gauge- invariant quantities ( for more 
details on how to determine the gauge dependence of perturbations in general, one can look up the review by Breuer 

§)■ 

When one chooses the I and n vectors of the unperturbed tetrad along the repeated principal null directions of the 
Weyl tensor, then 

00 = ^i = V># = i>£ = 

— = = = (7) 

By collecting the Newman-Penrose equations that relate 0oj'0i an d 02 with the spin coefficients and tetrad com- 
ponents of the stress-energy tensor, and linearizing about the perturbed values, Teukolsky gets (after some algebra) 
the following decoupled equation for the perturbed ipo : 

[{D-3e + e* - 4p - p*)(A - 4 7 + ^) 

- {5 + it* - a* - 3/3 - 4t)(<T + tt - 4a) - 3ip 2 ]itfj = 47rJ o (8) 

where 

T Q = (5 + 7T* - a* - 3/3 - 4t)[(Z> - 2e - 2p*)T l B m 
- (5 + it* - 2a* - 2P)T t f } + (D - 3e + e* - Ap - p*) 

x [{6 + 2tt* - 2/3)T£ - (D + 26 + 2e* - p* )T*J (9) 

Since the full set of NP equations remains invariant under the interchange I n, m «-» m* (this is the basis of 
the GHP method [Q), then by applying this transformation one can derive a similar equation for : 

[(A + 3 7 - 7 * + 4fi + n")(D + 4e - p) 

- (6* - t* + fi* + 3a + Att){5 - r + 4/3) - 3ip 2 ]i$ = 47rT 4 (10) 

where 

T 4 = (A + 3 7 - 7 * + 4/i + //)[(£* - 2t* + 2a)T r f m , 
- (A + 2 7 - 2 7 * + p*)T^ m ,] + (6* - T *+p*+3a + 4tt) 

x [(A + 2 7 + 2p*)T^ m , - (5* -r*+ 2(5* + 2a)T r f„] (11) 

In a similar way, one can define tetrad components of the electromagnetic field tensor 

$o = F^Pm", $i = \ F^n v + m^m"), $ 2 = F^m^n" (12) 

and get similar decoupled equations for <i>o and $2- One can try these ideas with neutrino and scalar fields also. 
So if one now writes the tetrads in Boyer-Lindquist |8) coordinates t, r, 9, <j) (after using the gauge freedom to set up 
the spin coefficient e = 0) || they become 

l» = [{r 2 + a 2 )/ A, 1, 0, a/A], n» = [r 2 + a 2 , -A, 0, a]/(2E), 

= [iasin9, 0, l,i/sin6]/(V2(r + iacosO)) (13) 
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where aM is the angular momentum of the black hole, £ 



i 2 cos 2 9, and A 



2Mr + a 2 (note that before 



one of the differential operators n^d^ was given the symbol A, but from now on it will have the conventional sense 
described here). With these expressions one can now write explicitly the spin coefficients and %p 2 - Then it turns out 
that one can write all the decoupled equations for test scalar fields (s = 0), a test neutrino field (s = a test 

electromagnetic field (s = ±1) or a gravitational perturbation (s = ±2) as a single master equation which is the 
famed Teukolsky equation : 



(r 2 
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a 2 sin 2 
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dr \ dr 
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d 2 j> 
d(j) 2 
i cosb 
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r — i a cosb 



dtp 
~dt 



[s 2 cot 2 e - s] ip = 4ttST 



(14) 



For the case that interest us, which is where the perturbations are to be interpreted as gravitational radiation that 
can be measured at infinity, the value for s = -2, ip = p~ 4 tpf , where p = — l/(r — iacosO) in the coordinates we are 
using, and T = 2p~^T^. 



As mentioned before, this equation turns out to be separable. If one writes the Teukolsky perturbative function as 
_ e -iujt e %m4> g^^ r ^ ftien the Teukolsky equation separates into a radial part and an angular part. The angular 



equation for S(9) has as a complete set of eigenfunctions the "spin weighted spheroidal harmonics" |J of weight s 
The radial part can be written as 



dr^ ' dr 



p 2 U)R = p 2 T 



(15) 



Inhere,p(r) = (r 2 -2Afr) _1 and the effective potential is given by U = (l-^)~ 1 [{Lur) 2 ~4iuj(r~3M)}~(l~l)(l+2) 1 
while T is the source term previously defined for s = —2. (This is for the simplified a = case. A slightly more 
complicated version depending on the value of a for the fully rotating case can be found as equation (4.9) of Teukolsky 

When the angular momentum parameter a goes to zero, the Kerr metric goes to the Schwarzschild one and the 
Teukolsky equation becomes the Bardeen-Press equation fTo| . 

There is of course considerable interest in computing the energy carried off by outgoing waves at infinity due to the 
evolution of the perturbations from some initial data. The non-trivial information about outgoing waves at infinity 
is carried by the ipf tetrad component. In principle, it is possible to use the solution for ipf to solve the complete 
Newman-Penrose set of equations for the perturbations in the metric. So, for outgoing waves with frequency u> 



-^(hfe-ihfj^ 



(16) 



Therefore, 



d 2 E^ _ lim r 2 ^ 
dtdfl ~ "°° I6tt 



B \21 



lim 



(17) 



III. CALCULATION OF PERTURBATIONS FROM INFALLING PARTICLES 



Both the Regge- Wheeler and the Teukolsky equation have been used extensively in one of the simple test cases for 
this perturbation formalism : that of a particle of mass p -C M falling into a stationary, isolated black hole. That 
leads later on to extensions like considering the deformation and internal dynamics of an infalling star, accretion disks 
p| , and hopefully the late stages of a black hole collision in a not-too-distant future. 

The first such calculation using a Green's function technique to integrate the recently derived Zerilli equation 
(the even-parity counterpart of the Regge- Wheeler equation) was done by Davis, Ruffini, Press and Price [O. They 
computed for the first time the amount of energy that was given out as gravitational waves by a particle falling radially 
from infinity into a Schwarzschild black hole, and they found that it radiated AE — 0.0104 p 2 /M in geometrized units. 
The radiation from the / =2 multipole dominates the spectrum and is peaked at ui — 0.32 A/ -1 , which just a little 
below the fundamental resonant frequency for the black hole. 
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Sometime later, Ruffini |l3|] treated the case of a particle falling radially from infinity but with non-zero initial 
velocity, with the result that the increase in radiated energy was minimal. More general treatments were attempted 
by Detweiler & Szedenits pj which examined infall trajectories with nonzero angular momentum. Considerable 
increases in the emitted gravitational radiation are seen as the normalized angular momentum of the trajectory 
J/nM increases from to close to AM (where the particle approaches a marginally bound, circular orbit). Increases 
in AE by a factor of 50 are found at the high J end. 

The first calculation involving a Kerr black hole was carried out by Sasaki and Nakamura Jlf| who considered 
a particle falling radially along the symmetry axis of the hole. Several studies (all in the frequency domain) have 
been carried out dealing with infall in the equatorial plane (and the effect of the rotational frame dragging) |L6| ], 
infall trajectories with finite angular momentum JL7| |, and quite a few detailed simulations of the gravitational waves 
emitted by a particle in a bound orbit around the black hole. O] 



IV. CONCLUSIONS AND CURRENT RESEARCH 



The Teukolsky equation is a powerful and very convenient way to deal with gauge invariant curvature perturbations 
of both the Kerr and Schwarzschild metrics. It has a very nice, separable mathematical structure and is amenable to 
robust numerical integration. Many interesting results have been derived leading up to the ideal situation of using it 
to treat the close limit of a generic rotating black hole collision and the evolution of gravitational waveforms from it. 

Up to now, basically all treatments have been based on the separability of the equation and calculate the energy 
and waveforms for the first few / multipoles of the spheroidal harmonics expansion once the radial part has been dealt 
with. However, for the purpose of detecting the gravitational waves from the inspiral collision of a binary black hole 
system using laser interferometers one would like to obtain the time integration of the full Teukolsky equation once 
we have started from reasonable initial data describing the two holes in close proximity to each other. 

Krivan et al. JL9| have devised a procedure to evolve perturbations in time from generic initial data using the 
Teukolsky equation. Their method analyzes the radiation at infinity by dealing with the s = —2 version of the 
equation. They avoid fully separating the Teukolsky function and use the ansatz 

ip = e' im V<I>(i,r*,(9) (18) 

Then the equation is rewritten as a first order matrix equation and numerically integrated. It has given encouraging 
results in treating scalar fields, scattering of gravitational waves and analysis of quasi-normal ringing and power law 
tails of the outgoing radiation. This demonstrated the feasibility of this numerical approach for the homogeneous 
Teukolsky equation. 

We are now working on the case of the infalling particle to again try to reproduce the Davis et al. results. The 
idea is to simulate the stress-energy tensor for a point particle 

!->(.) -/^^-.(r))* (19) 

by using very narrow gaussian distributions in place of the spatial Dirac delta function remaining after the integration. 
In here, z(t) is the geodesic trajectory of the falling particle of mass fj,, and u a is its four-velocity. Then one contracts 
with the null tetrad vector and inputs this into equation ([n]). Since this is a 2-d code that evolves the Teukolsky 
equation in an (r, 9) grid, one must carefully keep track of the <fi dependence and all its derivatives implicit in equation 
( |i"l| ) and then Fourier expand this piece of the source term to match the tf> dependence of the expression for i/j in 
equation (fl8|). 

We have had some encouraging preliminary results with the outlined procedure described here for the simplified 
case of a particle falling along a radial trajectory to a non-rotating black hole, and we are actively working to modify 
the code to tackle more general and interesting problems in perturbations of rotating black holes in the presence of 
infalling matter. 
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